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We address the long standing problem of recovering dynamical information from noisy acoustic 
emission signals arising from peeling of an adhesive tape subject to constant traction velocity. 
Using phase space reconstruction procedure we demonstrate the deterministic chaotic dynamics by 
establishing the existence of correlation dimension as also a positive Lyapunov exponent in a mid 
range of traction velocities. The results are explained on the basis of the model that also emphasizes 
the deterministic origin of acoustic emission by clarifying its connection to sticks-slip dynamics. 
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Adhesion continues to generate new directions of in- 
terest due to the wide ranging interdisciplinary issues in- 
volved and its technological importance. For instance, 
the recent surge in interest can be traced to its relevance 
to biological systems, in particular, the desire to design 
adhesive materials that mimic fibrillar adhesion inherent 
to biological species like gecko [1]. Despite the progress, 
day-to-day experience like acoustic emission (AE) during 
peeling of an adhesive tape has remained ill explained. 
This can be traced to the fact that most information is 
obtained from quasi-static or near steady state condi- 
tions and much less attention has been paid to nonequi- 
librium time dependent dissipative aspects of adhesion, 
and related phenomenon like friction (which is adhesion 
and wear) 0, ISj iM ^ ^1^° ^^1- kinetic and dynami- 
cal aspects involve interplay of internal relaxation time 
scales (determined by molecular mechanisms) with the 
applied time scale, they are important in a variety of 
situations that are subject to fluctuating forces such as 
flexible joints, composites, and even dynamics of cell ori- 
entation [5|. 

Dynamical information can be obtained using experi- 
ments on peeling of an adhesive tape mounted on a roller. 
These experiments show that peeling is jerky accompa- 
nied by a characteristic crackling noise p, The jerky 
nature is attributed to the switching of the peel process 
between two stable dissipative branches separated by an 
unstable one. ( The low and high velocity branches arise 
from viscous dissipation and brittle fracture respectively.) 
The negative force- velocity relation is common to many 
stick-slip situations, for example, sliding friction [1, Q 
and the Portevin-Le Chatelier (PLC) effect, a plastic in- 
stability observed in tensile deformation of dilute alloys 
Q, to name only two. In general, stick-slip dynamics 
results from a competition among inherent time scales 
[1)0, here, the viscoelastic time scale and the time scale 
of the pull speed. All stick-slip processes are examples of 
deterministic nonlinear dynamics. 

In contrast to stick-slip nature of peeling, the origin 
of AE ( even in the general context) is ill understood. 
Recently, we suggested that the energy dissipated in the 
form of AE can be modeled in terms of the local displace- 
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FIG. f : (a) A schematic of the experimental setup, (b) Plot 
of the scaled peel force function (^(w") as a function of . 



ment rate |10l |. A model relevant for the experimental 
set up that includes such a term reproduces major ex- 
perimental features of AE as also that of the peel front 
dynamics • The model also predicts spatio-temporal 
chaos for a specific set of parameters. Moreover, it is 
long believed that AE and stick-slip peel dynamics are 
related. But, establishing such a connection requires ex- 
tracting quantitative dynamical information from the AE 
signals which so far has not been possible largely due to 
the highly noisy nature of AE signals. Here, we show 
that deterministic dynamics governs the AE process by 
demonstrating the existence of chaotic dynamics using 
nonlinear time series analysis. The results are explained 
using a model that also provides insight into the connec- 
tion between AE signals and stick-slip dynamics. 

Retrieval of information about the underlying process 
is also important in the general context of AE as it is ob- 
served in a large number of systems like micro-fracturing 
process, volcanic activity [ll| , collective dislocation mo- 
tion [ij, [13I etc. However, most studies [ll|, [13], ex- 
cept Ref. [lj|, are simple statistical studies showing the 
power law distribution of AE signals as experimental re- 
alizations of self-organized criticality [14]. Even in Ref. 
[13 ] , the extracted fractal dynamics of dislocation gener- 
ated AE sources is aided by use of multiple transducers. 
However, the situation is more complex in peeling experi- 
ments as only a single transducer is used leading to scalar 
AE signals that are also substantially noisy making the 
intended task even more challenging. 

To verify the prediction of chaotic dynamics, we 
have performed peeling experiments of an adhesive tape 
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FIG. 3: (a) Correlation integral for pull velocity 4.8cm/s from 
d = 6 to 10. Dashed lines are guide to eye. (b) Lyapunov spec- 
trum of the AE signals for traction velocities V = 4.8cm/s. 



FIG. 2: (a) Raw and (b) cured AE signal for V = A.Scm/s. 
(c) Square of the amplitude ( in arbitrary units) for the data in 
(b). (d) model AE signal for V = 2.48 and m = 0.001 which 
is similar to (c) except for the magnitude of fluctuations. 



mounted on a roller driven at a constant traction velocity 
in the wide range 0.2 to 7.6 cm/s. A schematic of the ex- 
perimental setup is shown in Fig. [T^. An adhesive roller 
tape of radius R is mounted on an axis passing through 
O with a motor positioned at O' that provides a constant 
pull speed V. AE signals associated with stick-slip dy- 
namics are monitored using a high quality microphone. 
Signals were digitized at the standard audio sampling fre- 
quency of 44.1 kHz (having 6 kHz band width) with 16 
bit signals stored in raw binary files. For low pull speeds 
V, regular AE bursts are seen that correspond to stick- 
slip events separated by oscillatory decaying amplitude. 
With increasing pull velocity, the AE bursts become ir- 
regular and continuous as shown in Fig. [5^. There are 
38 data files each containing 1.2 x 10^ points. As in most 
experiments on AE, signals are noisy. 

Time series analysis (TSA) begins by unfolding the 
dynamics through phase space reconstruction of the at- 
tractor by embedding the time series in a higher di- 
mensional space using a suitable time delay [l5|. Let 
[x{k),k = 1,2,3, ••■,iV] be the AE signal with At as 
the sampling time. Then, d— dimensional vectors are de- 
fined by ~ [x{k), x{k + t), - ■ ■ , x(k + {d — 1)t)]; k — 
1 , • • • , [A^ — (d — l)r] . The delay time r is either obtained 
from the autocorrelation function or mutual information 
[l6j. Then, the chaotic nature of the attractor is quanti- 
fied by establishing the existence of correlation dimension 
and a positive Lyapunov exponent. 

The correlation integral defined as the fraction of pairs 
of points and whose distance is less than r, is given 
by C{r) = ^ E.,, - € - where e(. ■ •) is the 
step function and Np the number of vector pairs summed. 
A window is imposed to exclude temporally correlated 
points [l^. If the attractor is self-similar then, C(r) ^ 
r'^ , where v is the correlation dimension [l5| . Then, as d 
is increased, one expects to find convergence of the slope 
dlnC{r)/dlnr to a finite value in the limit of small r. 
In practice, the scaling regime is found at intermediate 



length scales due to the presence of noise. 

As the AE signals are noisy, we have used a modified 
Eckmann's algorithm suitable for noisy time series [l3]- 
Briefly, Eckmann's algorithm [18] relies on connecting the 
initial small difference vector ~ to evolved difference 
vector through a set of tangent matrices. The number 
of neighbors used is typically min[2d, d + 4] contained in 
shell size defined by inner and outer radii u and eg 
respectively. ( ti also acts as a noise filter.) The modifi- 
cation we effect is to allow more number of neighbors so 
that the noise statistics superposed on the signal is sam- 
pled properly. We impose additional constraints that the 
sum of the exponents be negative for a dissipative sys- 
tem, and also demand the existence of stable positive and 
zero exponents ( a necessary requirement for continuous 
time systems like AE) over a finite range of shell sizes ts ■ 
The algorithm works well for reasonably high levels of 
noise in model systems [17] as also for experimental time 
series (for details, see Ref. [13] )■ We have also repeated 
the analysis using the TISEAN package [l^ . 

The data sets are first cured using a noise reduction 
technique [iGj . Figs. [2^ and b show the raw and cured 
data respectively for V = A.Scm/s. Clearly, the domi- 
nant features ( the peaks shown by arrows) of the time 
series are retained except that the amplitude is reduced 
[l^ . Indeed, the two stage power law distribution for the 
amplitude of AE signals for the raw data are retained ex- 
cept that the exponent for small amplitudes is reduced 
(from 0.33 to 0.24) without altering that for the large 
amplitudes. The cured data are used to calculate the cor- 
relation dimension for all the data files. However, while 
raw data are adequate for calculating the Lyapunov spec- 
trum from our algorithm, cured data are required for the 
TISEAN package. To reduce the computational time, 
only one fifth of the total points are used. 

The autocorrelation time is 4 units in sampling time. 
A smaller value of r = 1 is used to calculate C{r). A log- 
log plot of C(r) for the pull velocity A.Scm/s is shown in 
Fig. [3] a for d = 6 to 10. A scaling regime of three orders 
of magnitude is seen with v ~ 2.65 ± 0.05. However, 
converged values of v ( using our method and TISEAN 
package) are seen, only for data sets for pull velocities 
from 3.8 to 6.2cm/s with ly in the range 2.6 to 2.85±0.05. 

Using our algorithm, the calculated Lyapunov spec- 
trum is shown in Fig. [3] b for V = A.Scm/s keeping 
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to = 0.065. Note that the second exponent is close to 
zero as expected of continuous flow systems. We have 
calculated Lyapunov spectrum for the full range of trac- 
tion velocities and we find (stable) positive and zero ex- 
ponents [19] only in the region 3.8 to 6.2cm/s, consistent 
with the range of converged values of v. 

As a cross-check, we have calculated the Kaplan- 
Yorke dimension D^y from the relation D^y = j + 

|A^^\|' ;SLi ^» > 0; J^ltl ^» < O. For the case shown 
in Fig. [lb, we find D^y = 2 + 1.5/1.6 = 2.94 consistent 
with v obtained from C{r) [l9|. Similar deviations are 
seen for other pull velocities. The Dky values obtained 
from the TISEAN package are uniformly closer to the 
values, typically Dky = i^-|-0.1. Finally, we note that the 
positive exponent decreases toward the end of the chaotic 
domain (6.2 cm/s). These results show unambiguously 
that the underlying dynamics responsible for AE during 
peeling is chaotic in a mid range of pull speeds. 

To understand the results, consider a recent model for 
peeling of an adhesive tape [l3|- In Fig. [T^, the distance 
OO' is denoted by / and the peeled length of the tape 
PO' by L. The angle between the tangent to the contact 
point P ( projection of the contact line PQ onto the plane 
of the paper) and PO' is denoted by 6 and the angle 
APOO' by a. From Fig. [TJi, we get L cos 6 — ~l sin a 
and L sin 9 = I cos a — R. As the peel point P moves 
with a local velocity v, the pull velocity is given hy V — 
V + ii + R cos 9 a. Defining u{y) to be the displacement 
with respect to the uniform 'stuck' peel front and defining 
v{y), 9{y) and a{y) at all points y along the contact line, 
the above equation generalizes to 



[V - v[y) - u{y) - R a{y) cob 9{y)\dy = Q. (1) 



where 6 is the width of the tape. As the contact line 
dynamics is controlled by the soft glue, we assume that 
the effective elastic constant kg along the contact line 
is much smaller than that of the tape material kt . This 
implies that the force along PO' equilibrates fast and the 
integrand in Eq. ([1]) can be assumed to vanish for all y. 

The basic idea of the model is that while stick-slip 
dynamics is controlled by the peel force function /(f), 
the associated AE is the energy dissipated during rapid 
movement of the peel front. We begin by defining dimen- 
sionless variable t = Uyt, with — (kt/bp) where p is 
the mass per unit width of the length L. Similarly, we de- 
fine u — Xd, I = l^d, L — L^d and R = R'^d using a basic 
length scale d = fmax/h, where f^ax = f{v„iax) is the 
maximum value of the peel force function f{v). We define 
the scaled peel force function by (j){v'^) = f {v^ {v)) / fmax 
(Fig. [1)3). Here, = v/vcCOyd and V" = V/vct^ud 
are the dimensionless peel and pull velocities respec- 
tively with Vc = v,nax/^ud- Using a scaled variable 
r — y/a, with a referring to a unit length along the 
peel front, the scaled kinetic energy can be written as 
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Here the first term represents the rotational kinetic en- 
ergy and the second term, the kinetic energy of stretched 
tape. C/ — {fmaxlktf{p/£) represents the relative 
strength of the two terms, where ^ is the moment of iner- 
tia per unit width of the roller tape. The potential energy 

= \S'J''X\r)dr+'-fST 
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with fco = {kgb'^ /kto?). The first term arises from the 
displacement of the peel front due to stretching of the 
tape and the second term due to inhomogeneity along 
the front. The total dissipation is the sum of dissi- 
pation arising from the peel force function 0(w*) and 
from the rapid movement of the peel front given by 

= \S'J''S<^{v^{r))dv^dr + \STl'^ 
spectively. 4)^{v^) is assumed to be derivable from a po- 
tential function $(«'') = j 4){v'^)dv^ . The second term 
denoted by Rae is the Rayleigh dissipation functional 
which is interpreted as the energy dissipated in the form 
of AE. The scaled 7„ is related to the unsealed dissipation 
coefficient r„ through 7„ = r„a;„/(fcta^). 
The scaled local form of Eq. [1] is 

X^{V' - v')vc + R'' \- {sin a) a. (2) 
Using Lagrange equations of motion, we obtain 
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dr. 



Equations ([2][S]) are solved using an adaptive step size 
stiff differential equations solver (MATLAB 'odel5s') 
with open boundary conditions. The nature of the dy- 
namics depends on the pull velocity , the dissipation 
coefficient ^y and Of. Cf depends on the roller inertia 
I = ^b (10^^ < I < 10^^) and the tape mass m — pb { 
0.001 < m < 0.1). 7„ ranges from 0.001 to 0.1. Other 
parameters are fixed at i?''=0.35, P =3.5, fco = 0.1 ( 
kt = lOOOiV/m) and N = 50. The (unsealed) peel force 
function f[v) preserves major experimental features like 
the values of fmax, and the velocity jump [g. 

The results reported are for m — 0.001 and 0.055, 
/ = 0.01 and low dissipation coefficient 7„ = 0.01. Phys- 
ically, low 7„, implies weak coupling between velocities 
on neighboring points on the peel front. Thus, local dy- 
namics dominates and hence more ruggedness leading to 
higher dissipation Rae (than for large 7„). Indeed, even 
for low y , the peel front breaks up into stuck and peeled 
segments (see Fig. [4^ for = 2.48 and also Ref. 10]). 
Hence, the acoustic energy dissipated Rae is noisy. 

Several qualitative features of the experimental AE sig- 
nals such as the change from burst to continuous type 
with pull velocity are displayed by Rae- The observed 
two stage power law distribution for the experimental 
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FIG. 4: (a)Stuck-peeled configuration for — 2.48 and m — 
0.001. (b)Tlie corresponding Lyapunov spectrum for Rae- 

AE signals is reproduced by the model. For instance, 
for the model signal in Fig. the exponent values are 
niE — 0.6 and 2.0 for small and large values respectively, 
consistent with the two exponents tua — 0.24 and 3.0 for 
Fig. (Note that energy Rae is the square of AE am- 
plitude.) Ref. [Tojl also reports a spatio-temporal chaotic 
state that corresponds to "edge of peeling picture" for 
high tape mass m = 0.1, / = 0.01 and low pull speeds. 
However, as experimental AE signals become chaotic as 
a function pull velocity (not studied in Ref. JUj), the 
correct quantity to analyze is the energy dissipated in 
the form of AE, RAsit) ( an average over the peel front). 

Following the embedding technique, we have analyzed 
the model AE signal RAE{t) and computed the correla- 
tion dimensions and Lyapunov spectrum for the entire 
instability domain. We find stable positive and zero ex- 
ponents for a range of values. A plot of the spectrum 
for m = 0.001 and V = 2.48 [to = 0.08) is shown in 
Fig. gh which gives Dky ^ 2 + 0.32/0.77 = 2.4 while 
ly — 2.2 ± 0.02 [l9|. Converged values of v ranging from 
2.2 to 2.7 {Dky in the range 2.4 to 3.0) are seen in the 
sub-interval 1.48 < 1^" < 6.48 of the instability along 
with stable positive exponents. Similar converged values 
of v for m = 0.055 ( ranging from 2.6 to 3.2 with Dky in 
the range 2.7 to 3.3) are seen in a mid range of V. The 
value of the positive exponent decreases for large V . 

Several conclusions emerge from the study. First, the 
presence of chaos in experimental AE signals supported 
by the model shows that deterministic dynamics is re- 
sponsible for AE during peeling. Second, the model also 
provides answers to questions raised by the TSA. For 
instance, the model shows that while stick-slip is con- 
trolled by the peel force function, acoustic emission Rae 



is controlled by the local kinetic energy bursts on the peel 
front generated during switching between the stuck and 
peeled states (Fig. H^). This mechanism provides insight 
into the transition from burst to continuous type of AE. 
At low pull velocities, V^, the number of stuck segments 
are few, each containing many spatial points (Fig. S^), 
with only a few large velocity bursts leading to burst 
type RAE{t)- With increasing the number of stuck 
segments increases (each containing fewer points) with 
a large number of small local velocity bursts that there- 
fore lead to continuous AE signals (similar to Figs. 3b, c 
of Ref.[l0|). Hence, the decreasing trend of the positive 
Lyapunov exponent with pull velocity observed in exper- 
imental signals can be attributed to peel front breaking 
up into large number of small stuck-peel segments. Thus, 
the model provides insight and clarifies the connection 
between stick-slip and the AE process. The work also 
addresses the general problem of extracting dynamical 
information from noisy AE signals. 

Our study has relevance to time dependent issues of 
adhesion, in particular, to failure of adhesive joints and 
composites that are subject to fiuctuating loads. Specif- 
ically, the analysis suggests that a larger value of the 
positive Lyapunov exponent (its inverse giving the time 
scale) implies higher dissipation and hence earlier failure. 
Thus, using acoustic emission technique to monitor AE 
signals in these cases coupled with the estimation of the 
largest Lyapunov exponent could prove to be useful. 

Many of these features are common to the PLC ef- 
fect. The effect attributed to pinning and unpinning of 
dislocations from solute atmosphere, is clearly, a distinct 
physical process from peeling. Yet, the negative force- 
velocity relation and the existence of chaotic dynamics 
in a mid range of drive rates are seen both in experi- 
ments and a model for the PLC eflFect as well d, H, [13] . 
Dynamically, the existence of chaotic dynamics as also 
the decreasing trend of the positive Lyapunov exponent, 
seen in both in the PLC effect and peeling, is the result 
of a reverse forward Hopf bifurcation (HB) (end of the 
instability) that follows the forward HB (onset) As 
chaotic window is seen in both cases, it is likely that it 
is a general feature in other stick-slip situations that are 
limited to a window of drive rates. 
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